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In order to guide the extrapolation of the mass of the rho meson calculated in 
lattice QCD with dynamical fermions, we study the contributions to its self- 
energy which vary most rapidly as the quark mass approaches zero; from the 
processes p — > ujtt and p — > irir. It turns out that in analysing the most recent 
data from CP-PACS it is crucial to estimate the self-energy from p — * irir 
using the same grid of discrete momenta as included implicitly in the lattice 
simulation. The correction associated with the continuum, infinite volume 
limit can then be found by calculating the corresponding integrals exactly. 
Our error analysis suggests that a factor of 10 improvement in statistics at 
the lowest quark mass for which data currently exists would allow one to 
determine the physical rho mass to within 5%. Finally, our analysis throws 
new light on a long-standing problem with the J-parameter. 
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I. INTRODUCTION 

As the lightest vector meson, the p is of fundamental importance in the task of deriving 
hadron properties from QCD. Within lattice QCD the ratio of 7r to p masses is often used 
as a measure of the approach to the chiral limit. For a long time lattice calculations were 
restricted to values of m^/rrip above 0.8. However, with the remarkable improvements in 
actions, algorithms and computing power, there are now lattice QCD results with dynamical 
fermions available for hadron masses with current quark masses such that m^/uip is entering 
the chiral regime. Nevertheless, in order to compare with the properties of physical hadrons 
it is still necessary to extrapolate the results to realistic quark masses 0. 

In the past few years there have been some very promising developments in our under- 
standing of how to extrapolate lattice data for hadron properties, such as mass 0, magnetic 
moments 0, charge radii || and the moments of structure functions J4|, to the physical 
region. In doing so it is vital to include the rapid variation at small pion masses associated 
with those pion loops which yield the leading and next-to-leading non-analytic behaviour. 
(This was crucial in arriving at a reasonable value for the sigma commutator ||, for ex- 
ample.) However, a formal expansion of hadron properties in terms of m n fails to converge 
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up to the region where lattice data exists. The crucial physics insight which renders an 
accurate chiral extrapolation possible is that the source of the pion field is a complex system 
of quarks and gluons, with a finite size characterised by a scale A. When the pion mass is 
greater than A, so that the Compton wavelength of the pion is smaller than the extended 
source, pion loops are suppressed as powers of m n /A and hadron properties are smooth, 
slowly varying functions of the quark mass. However, for pion Compton wavelengths bigger 
than the source (m^ < A) one sees rapid, non-linear variations. Phenomenologically this 
transition occurs at m n ~ 500 MeV, or m^/m p around 0.5 - the region now being addressed 
by lattice simulations with dynamical fermions. 

Another difficulty associated with the extrapolation of lattice results that needs further 
investigation is the discretisation of momenta inherent in all lattice calculations. In this 
regard we mention not only the finite lattice spacing but the fact that there is a minimum 
possible non-zero momentum available because of the finite volume of the lattice. This 
issue is absolutely critical to the interpretation of the recent CP-PACS data for dynamical 
fermions ||, in which a first result[] is reported at m^/nip ~ 0.4. As we explain in detail, 
the only reason that it is possible to measure the p mass there is that the calculation is done 
in a finite volume. We show that taking the finite lattice size and finite lattice spacing into 
account is a necessary requirement when extrapolating to the physical region. These effects 
become especially significant for the case of the p meson which has a p-wave, two-pion decay 
mode. 

In the next section we summarise the key physical ideas and the necessary formulas for 
extrapolating the mass of the p meson to the physical pion mass. This includes a discussion 
of the limiting behaviour at small and large quark mass. We then show the result of our 
fitting procedure and analyse the uncertainty in extracting the p mass at the physical point. 
We show that a factor of 10 increase in the number of gauge field configurations at the 
lowest quark mass presently accessible would be sufficient to determine the physical p mass 
to within 5%. In section [TTJ we discuss the consequences of this analysis for the J-parameter 
and conclude with a brief summary and outlook for the future. 

II. CHIRAL EXTRAPOLATION FORMULA 

The success of our earlier work concerning the extrapolation of the N and A masses |I| 
leads us to consider a similar approach to the latest two-flavour, dynamical QCD data on the 
p meson Once again our guiding principle is to retain those self-energy contributions 

which yield the most rapid variation with m n near the chiral limit - i.e. those terms which 
yield the leading non-analytic (LNA) behaviour and the dominant next-to-leading non- 
analytic (NLNA) behaviour. These processes are illustrated in Fig. [1]. The p — * ujtc term, 
shown in Fig. |l](b), yields the LNA contribution to the p mass. The p — > tttt term (Fig. 
|l](a)) not only yields the NLNA behaviour but, of course, the width of the p once m n goes 
below m p /2. 



Although CP-PACS finds no evidence of residual errors for the lowest mass point, they caution 
that it is premature to draw firm conclusions based on the present low statistics. 
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In order to evaluate these self-energy terms we take the usual interactions [@,[|: 



= J/jwnr P M • (?f X (9^7?) - (d^jf) X 7?) , (1) 

and 

C Ufm = g afm e, ua(3 (SPf) • tT . (2) 

These lead to the following expressions in the limit where the mass of the vector mesons (p 
and u, taken to be degenerate) is much bigger than the mass of the pion: 

y P _ fprnr f°° dk k 4 ul. w (k) . . 

™ 6tt 2 Jo w„{k){wl{k) -p2/ 4 ) ' {) 

yp _9lpJh f 00 dk k A ul w (k) 

12tt 2 Jo wl(k) ' 1 ' 

In analogy with the heavy baryon limit, we have neglected the kinetic energy of the heavy 
vector mesons. Here ££ w and correspond to the processes shown in Figs. [j](a), and 0(b), 
respectively. The pion energy is given by w n (k) = \Jk 2 -\- m 2 , and u nn and dipole 
form factors governed by a mass parameter reflecting the finite size of the pion source. In 
the chiral limit these have the standard LNA and NLNA behaviour (independent of the 
forms chosen for u nw and u nul ) : 

Pmr 4 

^-VttInLNA = _ 4 7r 2^2 m ^ M"^) , 
2 

^ttJlNA — 247J- ' ^ 

while they are suppressed as inverse powers of m n once m n is comparable with the dipole 
mass parameter .0 Finally, the p — > 7T7t term contains the unitarity cut for < p p /2 (as 
well as an imaginary piece determined by the width). 

The formal solution to the Dyson-Schwinger equation for the p propagator places the 
self-energy contributions in the denominator of the propagator and thereby modifies the p 
mass as 
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mp = ymg + £ 

~ m + o 6 

2m 

where £ = + ££ w and the bare mass, m , is taken to be analytic in the quark mass. 
Guided by the lattice data at large we will take m to be cq + 02m 2 . 



2 Note that all masses (e.g. the p mass, p, p ) and coupling constants should, in principle, be 
evaluated in the chiral limit. However, as the variations from the physical values are typically of 
the order 10%, we use the physical values. 



3 



71 CO 

fa) (b) 
FIG. 1. The most significant self-energy contributions to the p meson mass. 



The dipole form factors are denned as 




(7) 
(8) 



where p n and \l p are the physical masses of the 7r and p mesons. The normalisation of u n7T 
is chosen to be unity at the p pole and the coupling constant, f p7T7T = 6.028, is chosen to 
reproduce the width of the p (i.e., the imaginary part of the self-energy). In the p — > uj-k 
case we take g Ufm = 16 GeV -1 [pj] . The m% dependence of the self-energies of (|3|) and ([|) is 
shown in Fig. by the dot-dash and dashed curves respectively. The interesting behaviour 
of the p — > 7T7T self-energy has been noted in many earlier works. In the context of lattice 



QCD it was discussed by DeGrand [12] and by Leinweber and Cohen [10| and most recently 



by Szczepaniak and Swanson ]T3|]. Other studies have looked at the self-energy as a function 
of p 2 (invariant mass of the vector meson) for mixed [I4|-16|. 

Finally, the lattice data alone cannot separately determine A n7T and A nu . In order to 
constrain them we have therefore made the reasonable, physical assumption that the size 
of the source of the pion field should be the same regardless of whether the intermediate 
state involves an u or a 7r. Thus we require that is chosen so as to reproduce the same 
mean-square radius of the source as would be generated by the choice of A noJ . Equating the 
mean-square radii results in the following relationship: 



A, 



2JAI. 



(9) 



An alternative procedure, which could be imposed in future analyses would be to constrain 
the difference in the meson self-energy terms to reproduce the observed p — uj mass difference 

EM7J- 



A. Fitting Procedure 

As we hinted in the introduction, the fact that CP-PACS is able to extract a measurement 
of the p mass at m^/rrip < 0.5 is at first sight extremely surprising. Once the p is able to 
decay one would expect to measure not the p mass but the two-pion threshold. The origin 
of this result is the quantisation of the pion momentum on the lattice and in particular the 
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FIG. 2. Variation with pion mass of the self-energy contributions to the p meson, Eqs. (j3|) 
and d|]), for a dipole form factor with A nuj = 630 MeV . The solid points indicate the value of 
the self-energy when calculated at the discrete momenta allowed on the lattices considered in this 
investigation. The difference between the curves and points is an indication of the physics missing 
because of finite lattice size and spacing. 



fact that the lowest (non-zero) pion momentum available is 2n/aL where is L is the spatial 
dimension of the lattice. For the relatively small lattice used by CP-PACS at the lowest 
pion mass this corresponds to more than 400 MeV/c momentum. This is why the p remains 
stable. 

Motivated by Eq. (|B|), and wishing to preserve the correct leading non-analytic behaviour 
of the self-energies, we have chosen to fit the p mass with the simple phenomenological form: 

m p - c + c 2 m, + ————— . (10) 

Given the constraint relating A^ and A W7r , this involves three adjustable parameters. At 
large m v the self-energies are suppressed by inverse powers of m n and the p mass becomes 
a simple linear function of (or the quark mass). 

In the finite periodic volume, of the lattice, the available momenta, k, are discrete 



where L M is the number of lattice sites in the p direction, and the integer n M obeys 



Therefore to simulate the calculations that are done on the lattice, we should replace the 
continuous integrals over k in Eqs. (§) and @ with a discrete sum over \k\. However when 
| A; | is zero, the case of a pion emitted with zero momentum, the integrands vanish, and hence 
do not contribute to the self-energy. In fact there is no contribution to the self-energies until 
kfj, = ±27r/aL /i . There is therefore a momentum gap on the lattice for p-wave channels, 
produced by this discretisation of momenta. 

We have investigated this momentum dependence by evaluating the self-energy integrals, 
Eqs. (0) and (^), by summing the integrand at the allowed values of the lattice 3- momenta 

where the k^ are defined by Eqs. (|TT| ) and ( |12"D and V is the spatial volume of the lattice. The 
results for the self-energy contributions are presented in Fig. |2|. The self-energy calculated on 
the lattice (the solid circles and triangles) differs little from the full self-energy calculation 
in the high quark mass (ni^) region. Furthermore, the effect in the p — ► ujtt self-energy 
contribution is also small at low pion mass. The biggest change is in the p — > titt self-energy 
calculation, at lower quark mass. This is the region in which one might expect the biggest 
corrections because one is approximating a principal value integral on a finite mesh. This 
change in behaviour, particularly at the lowest data point (m 2 « 0.1 GeV 2 ), indicates that 
the 7T7T self-energy contribution is significantly understated in the lattice simulations. Upon 
calculating the full self-energy contribution via the continuous integrals, the magnitude of 
the self-energy is increased by about 10 MeV, which is 30% of the self-energy contribution 
at this point. These results for and E£ w are used in Eq. (p~0|) to fit the lattice data. 

Recent dynamical fermion lattice QCD results are presented in Fig. |3|. The scale param- 
eters relating the lattice QCD results to physical quantities have been adjusted U by 5% for 
the CP-PACS and UKQCD results. The effect is to increase the p mass from CP-PACS and 
decrease the mass from UKQCD providing better agreement between the two independent 
simulations. As the x 2 of the following fits is dominated by the CP-PACS data, we focus on 
this data set. 

Our fits using Eq. ([[(]) are based on the lowest five lattice masses given by CP-PACS. We 
selected the lowest lying masses because to move further away from the chiral limit would 
necessitate additional terms beyond the first two analytic terms of Eq. (0). The results of 
the fit are shown as the open squares in Figs. ||, |], and |. The parameters of the fit Cq, C2, 
and A nuJ , are then used in an exact evaluation of Eq. (|TCi| ) using the full integrals in Eqs. 
(§) and (f|). This result is illustrated by the solid lines in Figs. |] and |^. We note that the 
value A^ = 630 MeV for the best fit results in a softer form factor than one might expect. 
We do not consider this to be of significant concern in the present work because, as we shall 
discuss below, the current lattice results at low are not precise enough to constrain the 
chiral behaviour. 

It is interesting to note the similarity of the results to those of Ref. [r0|. There it 



was found that fitting quenched lattice data with a linear extrapolation, and improving the 
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FIG. 3. Vector meson (p) mass from CP-PACS § (filled circles) and UKQCD @ (open circles) 
as a function of m 2 . The dash-dot curve is the naive three parameter fit, Eq. (|i~3|). The open 
squares (which are barely distinguishable from the data) represent the fit of Eq. ((TTJ) to the data 
with the self-energy contributions calculated as a discrete sum of allowed lattice momenta. We have 
used a dipole form factor, with = 630 MeV. The solid curve is Eq. (]i~3| ) with the parameter C3 
fixed to the value given by chiral perturbation theory. 



extrapolation by adding on the p — > im effects, predicted essentially the same physical mass, 
but that the chiral behaviour was significantly different. 

For comparison we also show a popular three parameter fit, motivated by chiral pertur- 
bation theory: 



m. 



c + c 2 m^ + c 3 m m 



(13) 



This naive three parameter fit is illustrated by the dash-dot curve in Fig. |3|. However 
since the value of C3 in Eq. ( |l~3|) is commonly treated as a fitting parameter, we are not 
guaranteed that it has the correct value required by Chiral Perturbation Theory (%PT). 
The value for the best fit is found to be —0.21 GeV -2 . As outlined above, our expressions 
for the p self-energies have the correct LNA and NLNA coefficients by construction. Indeed, 
if the coefficient C3 is constrained to the correct valueQ (— g£ /487r = —1.70 GeV -2 ), the 
best fit possible by varying c\ and c% is shown as the solid line in Fig. [3]. As was also found 



3 In Ref. [O] the m n dependence of the LNA term to the p mass is given by 
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FIG. 4. Analysis of the lattice data for the vector meson (p) mass calculated by CP-PACS 
as a function of m 2 . The squares represent the fit of Eq. ( |l~0| ) to the data with the self-energy 
contributions calculated as a discrete sum of allowed lattice momenta. The solid curve is for 
continuous (integral) self-energy contributions to Eq. (13). We have used a dipole form factor, 
with optimal A nuj = 630 MeV. The shaded area is bounded below by a la error bar. The upper 
bound is limited by the constraint A^ > fi w as discussed in the text. 



in the case of the nucleon the lack of convergence of the formal expansion is such that 
it is not sufficient to fix the coefficient of the LNA term in a cubic fit to that predicted by 
xPT, as the resulting form will not fit the data. 

The importance of the accuracy of the lowest mass point cannot be overstated. We stress 
that CP-PACS emphasised the preliminary nature of the lowest data point, because of the 
relatively low statistics. Nevertheless, in order to prepare for future more accurate data, we 
have carried out a standard error analysis including this point and the results are presented 
in Fig. f|. The lower bound on the shaded area was found by increasing the minimum \ 2 
per degree of freedom of the fit by 1. We were unable to do this with the upper bound. The 
result is actually limited by the physics of the process. In the case of a dipole form factor 
this means A^ > fi n , and that is the upper limit we have shown here. 

It is not unreasonable to expect an improvement in the accuracy of the calculated lattice 



This would result in a value of the coefficient of 
value used here. 



T.71 GeV 2 , in excellent agreement with the 
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FIG. 5. The graph is as described in Fig. || except that the error bar on the lowest data point 
(m 2 w 0.1 GeV 2 ) has been reduced by a factor of \/T0. This equates to an improvement of 10 
times in the statistics, which we do not consider an unreasonable goal for the future. The dipole 
mass of the best fit is then = 660 MeV. The shaded area is bounded above and below by a la 
error bar. 



mass values, and as a Gedanken experiment we have explored the possibility of a ten-fold 
increase in the number of gauge configurations at the lowest pion mass. For the purposes 
of the simulation we did not change the value of the data point, but simply reduced the 
size of the error bar by vlO. As can be seen in Fig. [5] the improvement in the predictive 
power is dramatic. The uncertainty in the physical mass has been reduced to the 2% 
level. Additional improvement in the accuracy of the extrapolation would result from the 
availability of additional data in the low pion mass region. However, it must be noted 
that the provision of data around 0.2 GeV 2 and higher would probably not assist greatly in 
the determination of the dipole mass (A); it is primarily determined by points nearer the 
physical region. We present the parameters of these fits in Table |. 

We have examined the model dependence of our work by repeating the above fits with 
a monopole form factor. As can be seen in Fig. || the model dependence is at the level of 
15 MeV at the physical pion mass with current data, and at the few MeV level had the 
error bar been reduced by a factor of VTO. This reinforces the claim in Ref. [|l[] that this 
extrapolation method is not very sensitive to the form chosen for the ultra-violet cut-off. 
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FIG. 6. A magnification of the physical pion mass region of our extrapolation results. The 
solid and long dashed lines represent the best fit dipole and monopole results for a fit with the 
present accuracy of the lattice QCD results. The dash-dot and short dashed lines are the dipole 
and monopole results for a reduction in the error bar of the lowest lattice data by a factor of \/l0. 
The model dependence of the choice of form factor is 0(2%). 



III. J - PARAMETER 



A commonly perceived failure with quenched lattice QCD calculations of meson masses 
is the inability to correctly determine the J-parameter. This dimensionless parameter was 
proposed as a quantitative measure, independent of chiral extrapolations, thus making it an 
ideal lattice observable [ I9| . The form of the J-parameter is: 



J 



dm. 



m 



m K 



- m f 



=1.8 



tjik* ■ 



m 



K 



mi 



(14) 
(15) 



By using Eq. ([To] ) and the experimentally measured masses of the K (495.7 MeV), K* (892.1 
MeV), vr (138.0 MeV) and p (770.0 MeV) Lacock and Michael ||] determined 

J = 0.48(2) . 

However previous attempts by the lattice community to reproduce this value have been 
around 20% too small. In the case of quenched calculations this has been cited as evidence 
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FIG. 7. The solid curve is a plot of the value of the J-parameter as a function of m 2 obtained 
from Eq. ( |l4] ) and the best fit to the lattice results given by Eq. . The vertical dotted line shows 
the point at which the J-parameter is evaluated {mp/m^ = 1.8). The horizontal line displays the 
experimental value (0.48) plotted between the physical values of m 2 and rr? K . 



of a quenching error (see, for example the review in |2D|). It was noted by Lee and Leinweber 
| 2l| that the inclusion of the self-energy of the p- meson generated by two-pion intermediate 
states (excluded in the quenched calculations) acts to increase the J-parameter. 

In Fig. [F| we present the value of the J parameter obtained from Eq. ([14]) and our best 
fit to the lattice results using Eq. (0). The vertical dotted line indicates the value of m 2 
where the J parameter is to be evaluated, i.e. nip/m^ = 1.8. The horizontal dashed line, 
plotted between the values of the squares of the physical pion and kaon masses, shows 
the experimental estimate of the J parameter from flip]). This equation suggests that the 
evaluation of J may be approximated by the slope of the vector meson mass extrapolation 
between these points. The cusp shown in Fig. [7], associated with the cut in S^, suggests 
otherwise. We stress that while the detailed slope of the curve is parameter dependent, the 
presence of the cusp is a model independent consequence of the two pion cut in the rho 
spectral function. 

As a point of comparison we have also calculated J using the naive cubic chiral extrapo- 
lation, Eq. (0), described above. The results of our investigations are summarised in Table 
|. The value of the J parameter is similar for both fits as it is evaluated at m 2 ~ 0.22 GeV 2 . 
The effects introduced into the extrapolations by chiral physics do not begin playing a large 
role until m 2 falls below 0.2 GeV 2 . Had the J parameter been evaluated at m 2 = 0.19 GeV 2 
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Fit Form 


CO 


C2 


C3 




M p 


J 


2 

mi 


Cubic 


0.723 


0.668 


-0.207 




0.735 


0.44 (8) 


0.223(7) 


Dipolc 


0.776 


0.42732 




0.630 


0.731 


0.45 (7) 


0.225(4) 



TABLE I. Table of fit parameters Co, C2, C3, A^, the p-meson mass at fi^, the value of the 
J-parameter, and the pion mass at which the J parameter is calculated. All values are in appro- 
priate powers of GeV. The Cubic fit refers to Eq. (|l~3|) while the Dipole refers to Eq. ( |To| ) with 
a dipole form factor. We find that the error in the J-parameter is halved if the statistics on the 
lowest point are increased by a factor of 10. 



or 0.09 GeV 2 one would find perfect agreement with the linear ansatz of Eq. (|l~5"l). 

IV. CONCLUSION 

We have explored the quark mass dependence of the p meson including the constraints 
imposed by chiral symmetry. The pionic self-energy diagrams are unique in that they give 
rise to the leading (and next-to-leading) non-analytic behaviour and yield a rapid variation of 
the meson mass near the chiral limit. These are the lowest energy states with given quantum 
numbers that have significant couplings to the p-meson. Other meson intermediate states 
are suppressed by large mass terms in the denominators of the propagators, and also by 
smaller couplings. 

We find that the predictions of two-flavour, dynamical-fermion lattice QCD results are 
succinctly described by equation fllCf ) with terms defined in (|J) and (f|) for m n < 800 MeV. 
We have shown that our formula gives model independent results at the 2% level for the 
physical mass of the p meson. However, firm conclusions concerning agreement between the 
extrapolated lattice results and experiment cannot be made until the systematic errors in 
the extraction of the scale of masses can be reduced below the current level of 10% and 
accurate measurements are made at m n ~ 300 MeV or lower. 

We have also calculated the J parameter by directly evaluating the derivative of our mass 
extrapolation formula. We find that the empirical estimate based on differences of meson 
masses misses important non-analytic effects in the derivative of m p with respect to m 2 , as 
illustrated in Fig. [<]. 

Finally we have investigated the effects of an improvement in the statistics of the lattice 
data. Present lattice data is not yet sufficiently precise to independently constrain the 
behaviour near the chiral limit. With the best data available one finds a p-meson mass of 
731 MeV with la bounds at 675 and 1062 MeV. One could constrain the bounds by using 
phenomenological guidance for the form factors, but we would prefer to wait for better 
lattice data. Figure ^ suggests that the p-meson mass could be known to within 5% in the 
very near future. 
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ADDENDUM 



Since the submission of this manuscript the CP-PACS collaboration has released a 
preprint [25|, with work showing J as a function of mass. We note that their analysis 
does not address the chiral physics studied here. As a result, their curves will omit the 
general feature of a cusp in the J parameter as discussed in this manuscript. A similar com- 
ment applies to the MILC collaboration preprint [2^]. We look forward to seeing a similar 
analysis to that presented here applied to these new simulation results. 



APPENDIX: 

In this appendix we present the evaluation of the leading non-analytic terms of the E£ w 
and self-energy contributions to the p-meson mass. By the definition in Eq. all the 
non-analytic behaviour is contained in these two terms. 

We note that the form of the self-energy contribution from p — ► -kuj is the same as that 
for the process a^N discussed in Ref. [|I|. Using the results found in that paper we can write 
(for the choice of a sharp cutoff (9 (A — k)) for the form factor u nuJ ) 

= [rnl arctan (A) + ^ - Am^) . (Al) 

The chiral behaviour of this expression is obtained by expanding it in m n about m n = 
(the chiral limit). We find that in this limit 

v " ^ Ami + l< - \ml + 0(ml)) , (A2) 



12tt 2 \ 3 77 2 A 

3. 



with the leading non-analytic term being of order m 

2 

^ttJlna - 247T t ^ ' 

The p — > 7T7T self-energy contribution is slightly more complicated. If we again choose a 
^-function for the form factor we can analytically integrate Eq. (0) giving 

f 2 l 



6tt 2 2(/i p /2) 

'A - (p p /2) + y/A* + ml\ (A + (p p /2) + ^ A 2 + ml 
arctan | , | — arctan -j^=^=?= 

\ yjm% - M2f 
m+ (/Up/2) 



arctan 




X-(^/2)V \Jml-(p p /2) 2 



I A 2 + m 2 + A^ 

-(3m^-2( /Up /2) 2 )( /Up /2)ln( V -JL ) - A(^/2) ^ A 2 + m 2 ) , (A4) 
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where A regulates the cut off of the integral. The region in which we are interested (the 
chiral limit) has m n < (/i p /2). Thus the arguments of the arctans are complex. We use the 
relationship 

i f 1 — iz\ 

arctan(z) = - ml — I , (A5) 



2 \l + iz 

to rewrite this expression in terms of logarithms with real arguments. Collecting the 
logarithms together results in the following expression for the p — ► 7T7r self-energy, for 
77V < (/Up/2): 



s^p _ fpmr ' f // /n \2 2\ 3 / 2 



2 

on 

6n 2 2( 

m 2 K - (/V2) 2 ) + A\ml - 2(/, p /2) 2 ) - 2A(^/2)^/(A* + ml){{^/2) 2 - ml) 



In 



ml(A* + ml - (/i p /2) 2 ) 

, 0X ( Jk 2 + ml + A\ , ] 

-(3m 2 -2( /Jp /2) 2 )( /Jp /2)lnl V _^ j _ A(/x p /2)^/A* + m 2 . (A6) 

Looking at just the lowest order, non-analytic, terms in the expansion about m n = we 
have 

E -i- = J hw^) (( 2( ^ /2)3 - 3( ^ /2K + \ Si) 

+ (3m 2 - 2( /Jp /2) 2 ) (fi p /2)) \n(m n ) 



J OTV7T 



pTVTT 4 



4ttV 2 



m w ln(m*-) , (A7) 



which is the result given in Eq. (|5]). 
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